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Abstract 

A novel method of determining the mass of the W boson in the W + W~ — ► ivl'v' channel is 
presented and applied to 667 pb _1 of data recorded at center-of-mass energies in the range 183- 
207 GeV with the Opal detector at LEP. The measured energies of charged leptons and the 
results of a new procedure based on an approximate kinematic reconstruction of the events are 
combined to give: 

M w = 80.41 ± 0.41 ± 0.13 GeV, 

where the first error is statistical and the second is systematic. The systematic error is dominated 
by the uncertainty on the lepton energy, which is calibrated using data, and the parameterization 
of the variables used in the fitting, which is obtained using Monte Carlo events. Both of these 
are limited by statistics. 
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1 Introduction 



At LEP2 energies, four-fermion final states with two charged leptons and two neutrinos are dom- 
inated by the pair production of W bosons. The probability of both W bosons decaying into 
leptons is about 10%. These events are characterized by an acoplanar pair of charged leptons 
in the final state fij. Systematic uncertainties due to colour reconnection, Bose-Einstein corre- 
lations and hadronization modeling, which plague the W mass measurement in other channels, 
are absent. The full reconstruction of the W + W~ — ► Ivl'v' events is not possible due to the 
presence of at least two unobserved neutrinos in the final state. Therefore the W mass is deter- 
mined from the energy spectrum of the charged leptons and from an approximate reconstruction 
of the event referred to as the pseudo-mass method. 

This paper presents a measurement of the W mass M\y in the fully leptonic channel with the 
Opal detector, using an unbinned maximum likelihood fit. The data sample has an integrated 
luminosity of 667 pb _1 recorded at center-of-mass energies between 183 GeV and 207 GeV. 

A general description of the reconstruction of W + W _ — ► ivl'v' events can be found in 0. 
The end-points of the leptonic energy spectrum in W + W~ — ► ivl'v' events depend on the W 
mass. Neglecting the masses of the charged leptons and the finite width of the W boson, the 
energy of the charged leptons can be written in terms of M\v as []: 

where s is the square of the center-of-mass energy and B\ is the angle between the lepton 
direction measured in the W rest frame and the direction of the W in the laboratory frame. The 
latter is not known, so the W mass is determined mainly by the endpoints of the distribution, 
which correspond to cos#| = ±1. In practice, however, the end-points of the distribution are 
smeared considerably by the width of the W boson, by initial state radiation and by the detector 
resolution. These effects weaken the sensitivity of this variable. 

In order to include the information from the angle between the two leptons and the correlation 
between their energies, a second observable based on an approximate kinematic reconstruction 
of the event — the so-called pseudo-mass — is used in the analysis. The complete reconstruction 
of the event would require the determination of the four-momenta of both charged leptons and 
the two neutrinos, 12 quantities in total assuming the masses are known. The four-momenta 
of the two charged leptons are measured; therefore, the full reconstruction depends on the 
determination of the momenta of the two neutrinos. Assuming four-momentum conservation 
and equal masses for the two W bosons, five more constraints can be obtained. An additional 
arbitrary constraint is imposed to allow the "reconstruction" of the event. By assuming that 
both neutrinos are in the same plane as the charged leptons, the kinematics can be solved to 
yield a W pseudo-mass, which is quite sensitive to the true W mass. Due to a twofold ambiguity, 
two solutions are found for this variable: 

Ml = - —^({PVi?-Qvt)fat? + Vl) (2) 

±V(P? X Pl') 2 [(Pe> + P£) 2 (£beam " - (P + Q) 2 ]) , 

x The convention c=l is used throughout this paper. 



4 



where ^ 

-E'beam is the beam energy, E^, E^/ are the energies of the two charged leptons, mi, mg' are their 
masses and p^, ~pt' are their three-momenta. In the present analysis the charged leptons are 
considered to be massless particles because the values of the masses are very small in comparison 
with their energies. The sensitivity of both solutions was studied with Monte Carlo simulations, 
which showed that only the larger solution, M+, is sensitive to M\y. The distribution of this 
solution shows an edge where the value of the pseudo-mass is close to the true W mass. 

The pseudo-mass is used here for the first time in determining the W mass in leptonic W- 
pair decays. It is shown that the correlation between the mass determination from the single 
lepton energy spectrum and from the pseudo-mass is small, so that the mass measurement can 
be improved by combining the two methods. For the measurement presented in this paper 
the total error is dominated by the data statistical error. The largest contributions to the 
systematic uncertainty in this analysis are due to the uncertainties on the lepton energy which 
are determined with data and scale with the available statistics. In future experiments with 
larger statistics it should be possible to reduce both the statistical and systematic uncertainty 
significantly. 

2 Detector Description and Event Selection 

A detailed description of the Opal detector can be found in Q. The data sample used for this 
analysis corresponds to an accepted integrated luminosity of 667 pb™ 1 , evaluated using small 
angle Bhabha scattering events observed in the silicon tungsten forward calorimeter Q. 

The selection of the W + W~ — > Ivi'v 1 events and the identification of electrons and muons are 
described in For this analysis events with two charged leptons are required. In view of the 

experimental problems associated with reconstructing r-decays, the present analysis is restricted 
to events which contain at least one W decaying into an electron or a muon for the leptonic 
energy method (section 3.1.1) and events in which neither of the two W bosons decays into a 
tau for the pseudo-mass method (section 3.1.2). Table [l] summarizes the signal efficiencies and 
purities Q. The dominant background sources are Wez/, ZZ and dilepton events. These sources 
comprise 40%, 30% and 15% of the total background cross-section for the W + W~ — > IvHv' 
events @. A total of 1101 events for the whole range of center-of-mass energies are observed in 
the data. Table || summarizes the observed numbers of events for each center-of-mass energy 
together with the corresponding integrated luminosities. 

2.1 Monte Carlo Event Generators 

At each center-of-mass energy, samples of signal events were generated with five different W 
masses: 79.33 GeV, 79.83 GeV, 80.33 GeV, 80.83 GeV and 81.33 GeV. In the range of 183-189 
GeV, samples corresponding to the central value of 80.33 GeV were generated according to the 
CC03 ^ diagrams with KoralW ||, e + e~ — > ZZ and Zee background events were generated 
with Pythia [|l| and Wez; with KoralW |J . For the other W masses and for center-of-mass 
energies in the range of 192-207 GeV KoralW was used both for the signal and for these 

2 The leading order W + W~ production diagrams, i.e. the i-channel v e exchange and the s-channel Z°/7 ex- 
change, are referred to as CC03, following the notation of §1. 



four-fermion backgrounds. Other backgrounds were simulated by RADCOR [[Tl[] for multi-photon 
final states, KoralZ and Bhwide [ 13 1 for dilepton final states and two-photon events were 
simulated with Vermaseren Q and Herwig Q. Finally Yfsww3 Jl6| samples were used 
to perform systematic studies related to photon radiation in W + W~ events. 



2.2 Classification of Events 

The momentum resolutions obtainable with Opal for electrons and muons differ significantly 
at high energies. For electrons, the preferred measurement uses the electromagnetic calorimeter 
energy information which has a resolution of approximately 3% at 45 GeV, whereas for muons 
at the same energy, the charged particle momentum determined in the central tracker has a 
resolution of approximately 8% [jl7|j . Therefore, to maximize the sensitivity of the measurement, 
electrons and muons are treated separately by defining different classes of events for the leptonic 
energy and the pseudo-mass analyses, depending upon the identified lepton flavours. Identified 
taus are rejected since their energy cannot be determined. 

In the case of the leptonic energy the information of both charged leptons is used inde- 
pendently. For this variable two different classes are defined: the first class contains leptons 
identified as electrons and the second class contains leptons identified as muons. To increase the 
number of leptons used in the analysis, the higher energy identified electrons or muons in events 
selected as either ev e Tv T or ^u^ti> t are also used. 

For the pseudo-mass three classes of events are defined as follows: 

1. Events consisting of two leptons identified as e. 

2. Events consisting of two leptons identified as \i. 

3. Events consisting of two leptons identified as an e and a ji. 

Table |3] shows the classification of events for the leptonic energy and the pseudo-mass. The 
number of leptons selected for the leptonic energy method, and the number of events selected 
for the pseudo-mass method are summarized in Table |2[ 



3 Extraction of the W Mass 

3.1 Unbinned Maximum Likelihood Method 

The extraction of the W mass was performed by a simultaneous maximum likelihood fit to the 
leptonic energy and the pseudo-mass distributions using the data from center-of-mass energies 
from 183 GeV to 207 GeV. The fit is performed with the following product of likelihood 
functions: 

^=n(n^ x ii £ H (3) 

k=l \i=l 8=1 / 

where k denotes the eight center-of-mass energies included in the fit and i runs over the two 
classes defined for the leptonic energy (LE) and the three classes defined for the pseudo-mass 
(PM). 

The leptonic energy and the pseudo-mass distributions obtained from Monte Carlo simula- 
tions (including all background sources) are parameterized by appropriate analytic functions, 

f = f(P u ---,P N ;x), (4) 



fi 



which depend on N parameters, Pj, and on the variable x, which stands for either the leptonic 
energy or the pseudo-mass. For each parameter a linear dependence on M\v is assumed: 

Pi = + b] x M w (5) 

The coefficients and b\ are obtained by fitting the analytic functions to the leptonic energy 
and the pseudo-mass spectra generated at different W masses. The coefficient b\ is included only 
if it is not compatible with zero within one standard deviation (i.e. its numerical value is larger 
than the corresponding error); otherwise only a constant term is used. This parameterization is 
performed independently for each class defined for the sensitive variables. Details of the method 
are discussed in the following sections. 

3.1.1 Parameterization of the Leptonic Energy 

The leptonic energy spectrum is fitted with a function fy, which is the product of two Fermi 
functions and a linear function: 

h = fl x h x h (6) 

with 

h = —^ , /2 = P5(1 + P 6 XX), h = ~^ • (7) 

e p a + 1 e p 4 + 1 

Pi and P3 correspond to the points of inflection of the two Fermi functions, P2 and P4 to their 
widths, P5 is the constant term of the linear function and the slope is proportional to P%. 

The fit function depends therefore on six parameters. The overall normalization of the 
spectra can be used to eliminate P5. For each set of values of the Pi, P2, P3, P4 and Pq 
parameters, P5 is calculated so that the total function is normalized in the region taken to 
perform the fit. The lower limit is 10 GeV for all center-of-mass energies and the upper limit 
changes depending on the center-of-mass energy from 80 GeV for 183 GeV to 100 GeV for 
207 GeV. The regions at the end-points of the leptonic energy spectrum are sensitive to Mw- 
In terms of the parameterization these edges correspond to Pi and P3. These are the only 
parameters which are expected to change with the W mass while the other parameters, P2, 
P4 and Pq, are more sensitive to the W width and detector resolution effects. This has been 
confirmed by studying the linear dependence of each parameter on Mw- Figure [l] shows an 
example of a fit to the leptonic energy distribution. The spectrum is generated from a Monte 
Carlo sample with a W mass of 80.33 GeV at a center-of-mass energy of 189 GeV. 

A simultaneous fit of the parameters Pi(Mw), P2 and P3(M\y), P4 and Pq with the Monte 
Carlo samples for different values of Mw is performed. Parameters P2, P4 and Pq are assumed to 
be independent of Mw and a common value is determined for each of these. Pi and P3 describe 
the Mw dependence and a separate value is determined for each generated Mw- Figure [2] shows 
the fitted parameters Pi and P3 and the linear dependence on Mw for the case of electrons. 
These fits are performed independently for each class of events and each center-of-mass energy. 

3.1.2 Parameterization of the Pseudo-mass 

The analytic function chosen to fit the pseudo-mass spectra is the sum of a Fermi function and 
a constant function, with four free parameters: 



7 



f = Pl(^ E ^ + ^4). (8) 

\e p 3 +1 / 

Pi, P2 and P3 are proportional to the amplitude of the curve, to the point of inflection of the 
slope and to the width respectively and P4 is a constant term. As with the lepton energy, one 
parameter can be eliminated due to the constraint of the overall normalization. This procedure 
is used to eliminate the parameter Pi. For each set of values of the P2, P3, P4 parameters, Pi 
is determined to normalize the function within the fit range of 70 GeV to 90 GeV, which is the 
same for all center-of-mass energies. In terms of the above parameterization, the edge of the 
pseudo-mass distribution corresponds to Pj. This is the only parameter which is expected to 
show a clear dependence on M\v . Figure [3| shows an example of the function fitted to a simulated 
pseudo-mass distribution. Similar fits are performed for each class of events at all center-of-mass 
energies considered in the analysis. The edge around the generated W mass can be observed in 
the figure. The individual linear dependence of each parameter on M\v was studied in the same 
way as for the leptonic energy. In this case, a simultaneous fit of the parameters p2(Mw), P3 
and P4 with the Monte Carlo samples for the different values of M\v was performed. P3 and 
P4 are found to be independent of M\y; therefore a common value was determined for each. 
P2 describes the M\y dependence and an independent value was determined for each generated 
M\y- Figure || shows the fitted values of the P2 parameter and its linear dependence on M\v in 
the case of electron-electron events. 

3.2 Monte Carlo Studies 

The correlation between the fitted masses for the leptonic energy and for the pseudo-mass has 
been determined to be (11 ± 1)%, considering all center-of-mass energies together. This was 
evaluated by fitting the W mass separately for the leptonic energy and the pseudo-mass with 
an ensemble of 90 independent data-sized Monte Carlo subsamples generated with a W mass of 
80.33 GeV. 

The simultaneous fit method has been tested with 2000 data-size experiments derived by 
resampling from Monte Carlo simulations. The events of these subsamples were picked at ran- 
dom from the full Monte Carlo sample, whose signal part (generated at M\v = 80.33 GeV) 
corresponded to 90 times the data statistics. Multiple draws of the same events were allowed. 
Figure ||(a) shows the results of the resampling tests fitted with a Gaussian function. The pull 
distribution is shown in figure |5|(b). To compensate for the underestimation of the statistical 
error by the common fit due to the correlation between the leptonic energy and the pseudo-mass, 
the errors of the fit used in the pull distribution have been increased by a factor 1.11. 

3.2.1 Bias Test 

The method presented to extract the W mass is not expected to show any bias. Possible biases 
due to detector effects are eliminated since these are simulated with Monte Carlo methods. 
Biases related to the selection of the analytic functions which fit the leptonic energy and the 
pseudo-mass are negligible as long as the functions fit the variables properly. By using Monte 
Carlo samples generated with M\y = 79.33, 79.83, 80.33, 80.83 and 81.33 GeV, the bias and 
the linearity can be determined. First, the individual relations between the generated and the 
fitted masses are studied independently for the leptonic energy and for the pseudo-mass. The 
distributions show slopes and biases compatible with one and zero respectively, 0.958 ± 0.087 
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and -0.061 ± 0.054 GeV for the leptonic energy and 0.952 ± 0.086 and -0.031 ± 0.044 GeV for 
the pseudo-mass. Second, the relation between the generated and the fitted mass is studied for 
the simultaneous fit to both the leptonic energy and the pseudo-mass. This relation is shown in 
Figure || and found to be linear in a region of db 1 GeV from the central value of 80.33 GeV, with 
a slope of 0.983 db 0.063 and a bias value of —0.055 db 0.031 GeV. No corrections are included 
due to the slope or bias values since there is no evidence of a bias. 

4 Fit Results and Discussion 

The simultaneous fit to the data distribution of the leptonic energy and the pseudo-mass, com- 
bining center-of-mass energies from 183 GeV to 207 GeV, gives the following result: 

M w = 80.41 ± 0.41 GeV, 

where the quoted error is statistical only and has been scaled by a factor 1.11 to take into 
account the correlation between the two variables. The W mass values obtained separately are 
M w = 80.58 ± 0.52 GeV for the leptonic energy and M w = 80.20 ± 0.61 GeV for the pseudo- 
mass. Figure shows the comparison between the fit function and the data for the two classes 
defined for the leptonic energy at a center-of-mass energy of 207 GeV. The poorer resolution 
of the muon energy compared to the electron energy is visible in the fit functions. The regions 
around the end-points have sharper edges for electrons than for muons. Figure ^ shows analogous 
results for the three classes of events defined for the pseudo-mass at all center-of-mass energies. 
The greater sensitivity of electron-electron events relative to muon-muon or electron-muon events 
can be observed by comparing the fit functions around 80 GeV, i.e. in the region sensitive to 
Mw- 

5 Systematic Checks and Uncertainties 

The study of the systematic uncertainties on the W mass is described in this section and sum- 
marized in Table ||. The sources of systematic errors have been studied simultaneously for 
all center-of-mass energies. The total systematic error is calculated as a quadratic sum of the 
contributions listed below. 

5.1 Beam Energy 

The average LEP beam energy is currently known with a precision of about ± 20 MeV, varying 
slightly for different center-of-mass energies between 183 GeV and 207 GeV |i"q| . This leads 
to a systematic uncertainty of 11 MeV. The RMS spread in the LEP center-of-mass energy is 
around 240 MeV, varying slightly with the center-of-mass energy. The systematic error due to 
the uncertainty of the beam energy spread is 3 MeV. In addition the data at average center-of- 
mass energies of 205 GeV and 207 GeV were taken at a number of discrete energy points in the 
range of 204.6-206.0 GeV and 206.2-208.0 GeV, respectively. This leads to a bias of 12 MeV 
for which the fit result has been corrected. The average difference between the electron and the 
positron beam energies has negligible impact on the measurement of Mw- 
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5.2 QED Corrections 



The systematic error associated with uncertainties in the modeling of the QED corrections is 
estimated by reweighting events generated with Yfsww3, which has a more complete treatment 
of 0(a) QED corrections, according to KoralW probabilities. The corresponding difference of 
7 MeV is taken as systematic error. 



5.3 Detector and Resolution Effects 

The effects of the detector calibrations and the deficiencies in the Monte Carlo simulation of 
the detector response are investigated by varying the observed leptonic energy scales in the 
electromagnetic calorimeter (for the electrons) and in the central tracker (for the muons) over 



reasonable ranges [19|. The ranges used for the systematic variations depend on the polar angle 
and are determined from detailed comparisons of data and Monte Carlo utilizing both data 
recorded at 189-207 GeV and data collected at -,/s ~ Mz during 1998-2001. Due to differences 
observed in the resolution between the data and the Monte Carlo simulation during these years, 
correction factors for electron and muon energy are included in the reference samples. The 
difference between the measured and the generated energy is scaled by a factor 1.07 for the 
electromagnetic calorimeter and 1.06 for the central tracker. 

From the systematic studies of the lepton energy scale, the electromagnetic calorimeter 
energy scales and the central tracker scale are both known to 0.3%. These scale uncertainties 
dominate the detector related systematics in the W + W~ — ► Ivl'v 1 channel. For muons the scale 
uncertainty is 32 MeV. In the case of electrons, the scale errors are expected to be uncorrelated 
between the barrel and endcap regions and therefore the analysis for electrons was performed 
separately in both regions and combined in quadrature, resulting in an uncertainty of 85 MeV. 
The systematic uncertainty due to the non-linearity in the electromagnetic calorimeter and the 
central tracker is taken into account by varying the energy using a factor which changes linearly 
from —0.1% to 0.1% in the range 20-70 GeV. The corresponding systematic error is 43 MeV 
for electrons and 33 MeV for muons. 

The systematic error due to the uncertainties on the resolution is studied with Monte Carlo 
samples by comparing the fitted W mass with the nominal resolution to that with resolutions 
changed by 1.10 for electrons (with the information of the electromagnetic calorimeter) and 1.07 
for muons (with the information of the central tracker). The corresponding errors are 43 MeV 
and 12 MeV respectively. 



5.4 Background Treatment 

The background normalization is varied by ± 5%, which corresponds to the error in the cross 
section measurement for the leptonic channel @. The resulting change in the fitted M\v is 
11 MeV. 



5.5 Parameterization of the Sensitive Variables 

A total of 152 parameters are obtained for the leptonic energy and the pseudo-mass consider- 
ing center-of-mass energies from 183 GeV to 207 GeV for all the defined classes. From these, 
56 parameters depend on the mass of the W boson. To check the systematic error associated 
with the parameterization of the leptonic energy and the pseudo-mass distributions each pa- 
rameter is varied independently by ±lcr. This is repeated for all center-of-mass energies and all 
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classes defined for the leptonic energy and for the pseudo-mass. Adding the resulting changes 
in quadrature yields a total systematic uncertainty due to the parameterization of 51 MeV. 

Since the parameterization of the pseudo-mass does not describe a possible fall of the dis- 
tribution above 85 GeV we tried an alternative parameterization by multiplying the standard 
parameterization by a linear function. This changes the fit result by 4 MeV. 

5.6 Four-Fermion Effects 

Using events generated according to the CC03 diagrams instead of the full set of four-fermion 
diagrams results in a mass bias of 24 MeV. This has been determined by reweighting events 
produced with KoralW to the prediction of the CC03 matrix element. The fit result has been 
corrected for the mass bias due to the use of the CC03 Monte Carlo at center-of-mass of 183 GeV 
and 189 GeV. The remaining uncertainty due to the modeling of the four fermion final state 
is expected to be significantly smaller and has been neglected in the estimation of the total 
systematic uncertainty. 

6 Outlook for Future Experiments 

This analysis has demonstrated for the first time that the correlation between the two leptons 
can be used as additional information to the single lepton energy spectra in the determination 
of the W mass from leptonic W-pair decays. Due to the small correlation between the mass 
determination from the pseudo-mass and from the lepton energy the uncertainty is reduced 
significantly by combining the two measurements. The systematic uncertainty could be signif- 
icantly reduced in future high-statistics experiments. The uncertainties on the energy scale, 
resolution and linearity are determined with data and scale with the available statistics. The 
error due to the parameterization depends on the number of available Monte Carlo events and 
can be kept small by producing sufficiently large samples. 

The different sensitivities of the edges of lepton energy distributions and the pseudo-mass to 
the beam energy result in a reduced systematic error from the uncertainty in the beam energy in 
this analysis compared to the W + W~ — ► qq£v and W + W~ — ► qqqq channels |19(| . This will be 
important at future colliders, where the beam energy measurements will be one of the limiting 
systematic uncertainties. 

For the W + W~ — > qqlv and W + W~ — ► qqqq channels the hadronisation of the quarks gives 
rise to an important systematic uncertainty which can only be estimated by the comparison of 
different hadronisation models. The determination of the W mass in purely leptonic W decays 
is free from this uncertainty. As discussed, the total error in this channel should be competitive 
with the W + W~ — > qqiu and W + W~ — > qqqq channels with sufficiently high statistics. 

7 Summary 

The energy of the leptons produced in the reaction e + e~ — > W + W~ — > ivi'v' and the pseudo- 
mass, a variable which is based on an approximate kinematic reconstruction of the event, are used 
to measure the mass of the W boson. Both variables are combined in a simultaneous unbinned 
maximum likelihood fit, using data collected at center-of-mass energies between 183 GeV and 
207 GeV. The following result is obtained: 

M w = 80.41 ± 0.41 ± 0.13 GeV, 
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where the first uncertainty quoted is statistical and the second is systematic. The result obtained 
is consistent with previous measurements of M\y [19, 20jl- 
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Table 1: Signal efficiencies and purities for the leptonic events used in the analysis. 
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Integrated 
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Number of leptons 
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W+W~ -► IvHv' 


— leptonic energy — 


— pseudo-mass — 






events 






183 


57 


71 


77 


26 


189 


183 


278 


309 


80 


192 


29 


52 


51 


12 


196 


77 


144 


169 


32 


200 


74 


134 


144 


25 


202 


37 


82 


87 


16 


205 


82 


129 


141 


14 


207 
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Table 2: Integrated luminosities for the data from 183 to 207 GeV. The total number of events 
selected as W + W~ — > tvl'v' ', the number of leptons used for the leptonic energy and the number 
of events used for the pseudo-mass are shown. 
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Table 3: Classification of events used for the leptonic energy (LE) and the pseudo-mass (PM) 
analysis. 
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0.051 
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0.127 



Table 4: Summary of systematic uncertainties on the M\y measurement. 
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Figure 1: Fit to the leptonic energy distribution generated with M\v = 80.33 GeV for 
= 189 GeV. The crosses indicate Monte Carlo events and the shaded area shows the back- 
ground Monte Carlo. Only leptons tagged as electrons were used for this distribution. The fitted 
function has five free parameters and is a product of two Fermi functions and a linear function. 
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OPAL Monte Carlo 




Figure 2: Linear fit (Pi = b + b 1 x (M w - 80.33)) of the coefficients (a) Pi and (b) P 3 (points of 
inflection of the Fermi functions) at a center-of-mass energy of 189 GeV. The events chosen to 
perform this fit belong to the first class defined for the leptonic energy, which contains electrons 
only. Similar studies are performed for the second class and at all center-of-mass energies. 
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OPAL Monte Carlo 
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Figure 3: Fit to the pseudo-mass distribution generated with M\y = 80.33 GeV for 
-y/s = 189 GeV. The crosses indicate simulated Monte Carlo events and the shaded area shows 
the background Monte Carlo. Only events tagged as electron-electron were used. The fitted 
function has three free parameters and consists of a Fermi function plus a constant. 
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Figure 4: Linear fit (P2 = bo + bi x (M\y — 80.33)) of the coefficient P2 (point of inflection 
of the Fermi function) to M\v at a center-of-mass energy of 189 GeV. The events chosen to 
perform this fit belong to the first class defined for the pseudo-mass, which contains electrons 
only. Similar studies are performed for the second and third classes and at all center-of-mass 
energies. 
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OPAL Monte Carlo 
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Figure 5: (a) Fitted mass distribution from Monte Carlo for the simultaneous fit of the leptonic 
energy and the pseudo-mass at all center-of-mass energies. The test samples are generated at 
M\v = 80.33 GeV. (b) The corresponding pull distribution. The errors have been rescaled by a 
factor 1.11 to take into account the correlation between the leptonic energy and the pseudo-mass. 
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OPAL Monte Carlo 




Figure 6: Fitted versus generated W mass, showing the linearity of the unbinned method. The 
central value of 80.33 GeV is subtracted from all the masses. 



OPAL (207 GeV) 




Figure 7: Comparison between the data and the fit functions for the leptonic energy at a center- 
of-mass energy of 207 GeV. (a) electrons, (b) muons. 
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Figure 8: Comparison between the data and the fit function for the pseudo-mass at all center- 
of-mass energies, (a) electron-electron events, (b) muon-muon events, (c) electron-muon events. 
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